Spin-orbit interaction effect on transport of Dirac fermions in graphene 
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We study theoretically the quantum transport properties of the Dirac fermions with spin-orbit interactions 
(SOIs) in graphene by using the method of Schwinger proper time together with decomposition over Landau 
level poles and Kubo formula. The analytical expressions for both longitudinal and Hall conductivities are 
derived explicitly. It is found that, from some numerical examples, when the Rashba SOI is taken into account 
the Shubnikov-de Haas (SdH) oscillation peaks of the longitudinal conductivity versus the chemical potential 
are split, while the SdH oscillation of the longitudinal conductivity versus a external magnetic field exhibits a 
beating pattern. Furthermore, the Rashba SOI tends to suppress the quantum Hall effect in graphene. 
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I. INTRODUCTION 

Graphene has attracted a lot of attention because of its ap- 
pealing properties At low energy, owing to the specific 
band structure with the unique valley and neutrality separat- 
ing the hole states from the electron states^, graphene has led 
to the emergence of a paradigm of the relativistic condensed 
matter physics^'', where the relativistic quantum phenomena, 
some of which are unobservable in high energy physics, can 
be tested in the table-top experiment. The recent advances 
in fabrication technique have made it possible to produce 
graphite systems with a few layers or even a single mono- 
layer of graphene"" which propel graphene study into 
thriving high tide. 

The relativistic feature of the graphene represents the 
anomalous physical phenomena, such as anomalous quantized 
Hall effect, absence of the weak localization and existence of 
the minimal conductivity^. In addition to dissipative transport 
also supercurrent transport has already been observed"*. It is 
suggested that the graphene is a promising candidate for the 
spintronics and related applications due to the SOIs'"*'^. In 
a recent paper, Kane and Mele have studied the effect of the 
SOIs in graphene"''^, and find that the spin Hall conductivity 
is quantized in the absence of a magnetic field because of a 
gap produced by the SOIs. Depending on the relative strength 
of intrinsic and Rashba SOIs, it is further shown that the spin 
Hall conductivity can be zero or nonzero. Kane and Mele have 
roughly estimated the SOI scale. Following this routine, Yao 
et al}^ and Min et al. have found that this estimation is too 
big, and they have given some explicit expressions of the SOIs 
in graphene. 

In this work, we investigate the transport of Dirac fermions 
in graphene. The main purpose is focused on the effect of 
the SOIs on the transport on the basis of estimation of Yao 
et al. and Min et al. for the SOIs in graphene. Using the 
Schwinger proper-time method^S, decomposition over Lan- 
dau level poIes2iii and Kubo formula^^, we obtain some ana- 
lytical expressions for both longitudinal and Hall conductivi- 
ties. It is found that when the Rashba SOI is considered, the 
longitudinal conductivity as a function of the chemical po- 



tential deviates the linear relation at zero magnetic field. For 
nonzero magnetic field, the SdH oscillations are observed and 
the oscillation peaks in the longitudinal conductivity versus 
the chemical potential are split when the Rashba SOI is ap- 
plied, while the oscillation in the longitudinal conductivity 
versus the magnetic field exhibits a beating pattern. It is also 
shown that the Rashba SOI tends to suppress the quantum Hall 
effect in graphene. 

The rest of this paper is outlined as follows: In Sec. II, The 
model of a single layer of graphite (graphene) with the SOIs 
is established. In Sec. Ill, we derive the analytical expressions 
for both the longitudinal conductivity and the Hall conductiv- 
ity including the limits of these expressions at zero field. In 
Sec. IV, the corresponding numerical results and discussions 
are given. In Sec. V, a summery is presented. Finally, Some 
tedious algebra is included in the Appendix. 



II. MODEL FORMALISM 

The graphene is a flat monolayer of carbon atoms tightly 
packed into a honeymoon lattice. At low energy, it can be 
described by 2-1-1 dimensional relativistic field theory model. 
When the SOIs are included, the Lagrangian density of the 
system is given by 



(1) 



where *F - (^'ic,^k') is the eight-component Dirac spinors 
with ^K(K') = (*I'aT' 'J'bT' which describes the spin- 
related Bloch states residing on the atoms of the A, B sub- 
lattice at momentum K(K'), D - y^{d^ - ieA^) with y'^ 
(ju = 0, 1, 2) being Ax Ay matrices^^, e is the electron charge, 
Vf is the Fermi velocity, the external magnetic field B = V x A 
is applied perpendicular to the x - y plane and the coiTe- 
sponding vector potential is taken in the symmetric gauge 
A = {-By 12, Bx/2). In Eq. O, H, describes the SOIs that 
readiiL^ 



= Asod - y°s,) + Adiy's, - yVs,), 



(2) 
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where Aso is the intrinsic SOI parameter, Ar is the Rashba 
SOI parameter, and s is the spin variable. For Z? = 0, the 
corresponding energy spectrum are given by 



(3) 



For Aso > > 0, the system includes an energy gap of 
2(/lso - ^^r)- For < Aso < ^r, the energy gap closes. 

The Green's function of Dirac fermions described by the 
Lagrangian ^ in an external magnetic field can be expressed 
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+ D2 - i[D, Hs] 
Using the Schwinger proper time approaches, we obtain 



G(x,y) = expiie I A^c/z )G(x - y). 



ds 
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where f^v = ^/jA^ - dyA^,, = g'"' + [1 - eBs cot(eBs)] with cr'^" = 5(r''r''-r''r'')- It is clear that the symmetric gauge 



withg^" = diag(l, -1,-1), and 

B^ = i2Aso(^''s, + 2AR{(^'sy - S^oj' s, - S^xj^s,), (7) 

A = 2(4o + - /^J + 4^50/l«(/r' - /r' ^.v) (8) 



sets the factor exp(/e J^' A^iiiz'')=l in Eq. (|5]l. Therefore, Eq. 
Q becomes the functions of the difference x - y only. Using 
the expansions of the exponential function e^, one can show 
that 



r'<'°'*'5o+24)'[cos(^i) + //i- sin(^.s)]{cos(eBi) + y'y^ sin(eBi) + -(1 - /i,)[cos(^i) - cos(eBi)] 



.12iso^«sin(^i) 1 1 X 1,1 0^1 2,^^ w d mi 

I {ly'sy - yts^) + -(1 - r iz)r r [y sm(^i) - sm(eBi)]), 
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where ( = - 6^, ^ = ^j {2A As o^^r)'^ + ieB)^. Sub- 

stituting Eq. (??) into Eq. (|5]), applying Fourier transform 
in the Matsubara representation and using the decomposition 
method over Landau level poles^' we can derive the expres- 
sion 

G{icLi„, k) = Gi(/(y„„ k) + G2(/w,„, k) + Gi(iu),„, k) (10) 
after straightforward but somewhat complicated calculations. 
I 



The expressions of G,('Wm,k) (/=1,2,3) are very complicated 
and will be given in the appendix A. Whence, we can further 
obtain the retarded and advanced Green's functions by the an- 
alytic continuation G**'(w + /O, k) - G{iu),n co + /O, k) and 
G^^Ho) - iO, k) = G{i(jj,„ — > <jj - iO, k). When considering the 
influence of impurities, it is assumed that the scattering rate 
r on impurity is described phenomenologically by a constant, 
and then the Green's functions acquire the form 



G(*''^>(w,k) 



Gf^\oj±iTM) + Gf-^\a>-. 



iT, k)-HGf'*'(w±/r,k). 



(11) 
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In general, the scattering rate F, which is defined by T((jS) - 
-Iml/^{cL>), is a frequency-dependent quantity. It needs to be 
determined self-consistently from the Schwinger-Dyson equa- 
tions. The exact form of this equation actually depends on 



the impurity scattering fashion, such as short- or long-range 
scatterers. This kind of consideration have been made for 
graphene in Ref.-^. But in this paper, we mainly focus on 
the SOI eff'ect on transport, and neglect the exact form of in- 
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teractions between impurities and electrons. 



tained by analytical continuation of the Matsubara function 



III. ELECTRONIC CONDUCTIVITY 

The Kubo formula concerning the frequency-dependent 
electrical conductivity as a linear response function to an ex- 
ternal field can be written as^^ 



CTijiD.) 



Imn*.(Q + iO) 
Q ' 
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where i,j are the component indexes of coordinates, and 
llfj(oj) is the retarded current-current correlation function ob- 



1 

Uijiiojn) = - dTe'""\TrJi(T)Jjm, aj„ = InnT, (13) 
y Jo 



where V is the volume of the system, /3 is the inverse temper- 
ature, and Jiij) - J d^rji(T,r) with j,- - -ev/r*P(T, r)7'*P(T, r). 
Neglecting the impurity vertex corrections, the calculation of 
the conductivity reduces to evaluation of the bubble diagram. 
Then Eq. (fT2] l can be rewritten as 



nji^) = T^Re doj —rtr{[nF(cu)-nF(cj + n)][yG\aj + n,k)yjG^(aj,k)] 
ZttLI J {Zny 

+ npicj + Q)y'G^{u) + Q, k)yG'*(w, k) - np{u)y'G'^{cj + Q, k)y^G'^{(x>, k)]), 
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(14) 



where nf (cj) is the Fermi distribution function. Substituting 
Eq. ( fTTT i into Eq. (14), we can obtain the longitudinal conduc- 
tivity 



cr^„ - — T^Re duj- 

■^^ 71- j-co 4 
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and the Hall conductivity 
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-Ah(w), (16) 



where all the quantities on the right-hand side are calculated in 
the Appendix B. Eqs. ( fTSl l and ( fTSI l establish the fundamental 



basis for investigating the SOI eff'ect on the quantum transport 
properties of the Dirac fermions in graphene. 

In the limit of zero field, the Hall conductivity becomes 
zero. While for the longitudinal conductivity, using the 
asymptotic expansions 



iA(z) = ln^- 



l _ 1 1 1 

Yz " 12? 120z4 ^ V 



we arrive at 
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for Aso - 0. When T — > 0, >> F, the longitudinal 
conductivity can be further expressed as 



4e2 



lju|/r 



n l+(34/Fp)2 



(19) 



IV. RESULTS AND DISCUSSION 

To investigate numerically the behavior of electrical con- 
ductivity, we need to restore the whole model parameters in 



Eqs. ( fTSl l and (fTSI l. Thus, one should carry out the replace- 
ments: T — > ksT, eB — > tieBv^p. In the following, we mainly 
discuss the Rashba SOI effect on the transport properties since 
the intrinsic SOI is very small, while the Rashba SOI can be 
tunable by a perpendicular electric field. In Figs. 1 and 2, we 
show the chemical potential yu dependence of the longitudinal 
conductivity for the different Rashba spin orbit parameter Ar 
at zero or nonzero field. For zero field (see Fig. 1), one can 
see that when Ar - 0, the conductivity is proportional to |yu| 
and tends to the known quantum-limited minimal value 4e^//i 
at zero chemical potential'. For Ar + 0, there exists a thresh- 
old chemical potential [ic which increases with increasing Ar. 
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FIG. 1: The longitudinal conductivity cr„ measured in 2e^ jh units as 
a function of the chemical potential /j for the different values of Ar. 
We take S = OT, T = 3K, T = 5K, and Aso = O.OOIK. 



FIG. 3: The magnetic field dependence of the longitudinal con- 
ductivity measured in 2e- Ih units for the different values of Ar at 
= -600K. The other parameters are taken as Fig. [T] 




FIG. 2: The longitudinal conductivity measured in 2e^ jh units as a 
function of the chemical potential yu for the different values of Ar at 
B = IT. The other parameters are taken the same as Fig.[T] 



When the chemical potential is smaller than jic, the longitudi- 
nal conductivity becomes almost independent of yu; while for 
fi > He, the cr^v-A' curves recover the linear relation. This ten- 
dency agrees with Eq. ( fT9] l. For nonzero field case in Fig. 
2, we observe SdH oscillations of the conductivity due to the 
Landau-level crossing of the Fermi level^^. From Fig. 2, It is 
clearly seen that when Ar + 0, each oscillation peak is split 
into two implicit peaks, and the splitting peaks shift by Xr. 
This is due to the spin-orbit splitting of the Landau levels. 

The longitudinal conductivity as a function of the magnetic 
field B for the different Ar is shown in Fig. 3. For Ar - 0, 
the longitudinal conductivity decreases and intervals between 
the neighboring SdH oscillation peaks become large with in- 
creasing B, which reflects the fact that in the presence of the 
magnetic field only the transitions between neighboring Lan- 
dau levels contribute to electrical conductivity, while a further 



increase of the magnetic field leads to increasing of the dis- 
tance between neighboring Landau level, thus suppresses the 
transitions between them. When B is large enough, the con- 
ductivity becomes independent of B since the lowest Landau 
level is filled which is always below the Fermi level. These 
observations are quite consistent with the previous studies^^. 
In particular, it is interesting to note that when the Rashba SOI 
presents, the longitudinal conductivity exhibits the character- 
istic feature that the SdH oscillations are enhanced largely at 
certain positions, however damped at other positions. Such 
SdH as a beating pattern have been observed in two dimen- 
sional electron gas^^. From Fig. 3, one can find that the en- 
hanced positions and amplitudes of the SdH oscillations can 
be tuned by the Rashba SOI due to shift of one set of Landau 
level by Ar. 
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FIG. 4: The longitudinal conductivity measured in 2e- jh units as a 
function of Ar for the different magnetic field S at /j = -300K. The 
other parameters are taken the same as Fig. [T] 



5 



Figure 4 shows the longitudinal conductivity versus the 
Rashba SOI parameter /}« for the different magnetic field B. 
It is found that when the magnetic field is applied, the longi- 
tudinal conductivity as a function of Ar behaves as the oscilla- 
tion. It is because the Rashba SOI leads to the shift of landau 
level, the longitudinal conductivity shows a maximum each 
time a Landau level passes through the Fermi level of system, 
and a minimum when the Fermi level is situated between two 
Landau levels. 
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FIG. 5: The Hall conductivity measured in 2e^/h units as a function 
of the chemical potential yu for the different values of /l^ at S = IT. 
The other parameters are taken the same as Fig. [T] 

Figure 5 shows Hall conductivity as a function of the chem- 
ical potential fi for the different Ar. When Ar = 0, the Hall 
conductivity has a steplike structure as a function of fj., which 
reflect the quantum Hall effect. While the Rashba SOI opens, 
the Hall steps become naiTow and the step near = is split 
into two steps. It is observed that the Hall conductivity dis- 
plays peaks instead of a plateau at larger Ar. There is no Hall 
plateau in the cases of sufficiently strong Rashba SOI. This re- 
sult suggests that the Rashba SOI tends to suppress the quan- 

I 



tum Hall effect in graphene. 



V. SUMMARY 

We have investigated the effect of the SOIs on transport of 
Dirac fermions in graphene on the basis of amplitude estima- 
tion of Yao et al. and Min et al. for the SOIs. Using the 
Schwinger proper-time method, decomposition over Landau 
level poles and Kubo formula, we obtain the analytical ex- 
pressions for both longitudinal and Hall conductivities. It has 
been found that when the Rashba SOI is applied, the longitu- 
dinal conductivity versus the chemical potential deviates the 
linear relation at zero magnetic field. For nonzero magnetic 
field, the SdH oscillation in the longitudinal conductivity is 
observed, and each SdH oscillation peak is split into two peaks 
as the Rashba SOI is applied. While the oscillation in the lon- 
gitudinal conductivity as a function of the magnetic field ex- 
hibits a beating pattern with the Rashba SOI turned on. It is 
also shown that the Rashba SOI tends to suppress the quantum 
Hall effect in graphene. 
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APPENDIX A 

The Green functions Cf''^^ in Eq. (1 1) are given by 



- Af^\<.)] + (r% + r'k2)[Bf^U - Af '^V)], ^^^^ 

Gf ^V) = (Hs + 7°((o ± /r))7»iJAf '^'(w) - iyysgn(eB)Af^\oj)] - i(y'k2 - y^ki)y°s,sgn(eB)[Bf^\oj) 

- Af'\co)] + (y% + y%)y°sdBf^\co) - Af^\co)], ^ ^^^^ 

+ y k '(oj) + iy y sgn(eB)J\ '(a)))(iy Sy-y y s^) + {y ki + y k )(/g '{oj) 

+ iy'y^sgn(eB)jf-^\co))(iy'sy-yys.,)y's„ (A3) 



where Ai = /, + + /g - I-,, A2^h + Ulx + 4 - hlx^ A3 ^ h - 13 + I5 + h, A4 ^ h ~ hlx + h + hlx< 
Ji+ Jilx + Ja- hlx^ B^^ h- J2+h + J5, 84 = Ji- J^lx + Ja + Jelx^ in which 



j(R,A) 



(1 i\ = Y L„(2a)-L„-i(2a) , L„(2a)-L„_i(2a) -. 

n=0 (w±(T)2+104p+24-f+2n|eB| («±(r)2+104o+2^^+f+2n|eB| J' 

00 

r(^,A). _ 1 -a y. / lyr L„(2g)+L„^i(2a) L„(2a)+L„_i(2g) -, 

■'2 1^ ^ ^> L(a)±(T)2+104o+24-f+2«|eB| («±jT)2+104o+24+f+2n|eB|-'' 



2 ' '-(w±iT)2 + 104o+24-f-f+(2«+l)|eB[ (w±iT)2+104o+24+f+f+(2n+l)|eB| 

, ^(2«) , /.„(2a) 1 

(«±jT)2 + 104o+24-f+f+(2n+l)|eB| («±jr)2+104o+24+f-f+(2n+l)|eB|-'' 



r(fi,A). N _ l„-ay/_i\Kr ^(2a) L„ila) 

■'4 (.t^; ^\ ^) L(„±jT)2 + 104o+24-f-f+(2n+l)|eB| (w±iT)2+104o+24+f+f+(2n+l)|eB| 



(«±jT)2 + 104o+24-f+f+(2n+l)|eB| (w±(T)2+104o+24+f-f+(2«+l)|eB|-'' 

r(^,A)(- N _ 1 -a y / i -.Br Z-„(2a)-L„_i(2a) L^Clay-Ln^ijla) -, 

(^oi; 2^ ^^^'^ '-(w±!T)2+104p+24-f+2n|eB| («±;T)2+104o+24+f+2n|eB|-'' 

A _ 1 -a y / i yir L„(2g)+L„_i (2a) L„(2a)+L„_i (2a) -, 

'6 2*^ ^ ^> '-(w±iD2 + 104o+24-f+2n|eB| {oj±irf+m\lg+2\l+(+2n\eB\^' 

r(fi,A)/ A _ l„-ffy(-_l\nr ^„(2g) L„(2a) 

\a)) 2^ ^\ 1) '■(w±iT)2 + 104o+24-f-f+(2n+l)|eB| (w±iT)2 + 104o+24+f+f+(2n+l)|eB| 

, " ^(2«) /.„(2a) 1 

(«±jT)2 + 104o+24-f+f+(2n+l)|eB| («±jr)2+104o+24+f-f+(2n+l)|eB|-'' 

7(K,A). . _ i^-ayc_iYi| A,(2ff) Z.„(2a) 

(UJ) r^K i; L(^j_;pj2 + i(y2^^_|_2^2_^_^^(2«+l)|eB| (w±iT)- + l(W-5„+24+f+^+(2n+l)|eB| 



L„(2a) L„(2a) 



(w±ir)2 + 104o+24-f+f+(2n+l)|eB| (w±iT)2 + 104o+24+f-f+(2K+l)|eB| 

r(R4)/ ^ 1 -ay. 4L^_i(2a)+L„(2a)-L„_i(2a) 4L',_, (2a)+L„(2a)-Z.„_i (2a) 

7j 2^ ^Zj^"- '■(w±iT)2 + 104o+24-f+2nleBl (w±ir)2 + 104o+24+f+2n|eB| J' 



^(^-^V^)-) = l^-gy,- 1-)«r Z.,',(2a)+L',_,(2a) L;(2a)+£;,_,(2a) 

•'2 2*^ ■^^'^ '■(M±!T)2 + l(M-5o+24-f-f+(2n+l)|eB| (M±/r)2 + 104,j+24+f+f+(2n+l)|cB| 

Li(2a)+£;_,(2ff) L,',(2a)+L;_,(2a) , 



(w±ir)2 + 104o+24-f+f+(2n+l)|e£| (w±ir)2 + 104o+24+f-f+(2n+l)|eB|-" 
,(«^), . 1 -av^- IVT Z.i(2a)+L,;_,(2a) Lj(2a)+L;_,(2a) 

A (,^^; = 7^ L 



(w±iT)- + 104o+24-f-f+(2n+l)kB| (w±ir)2 + 104o+24+f+?+(2n+l)kB| 
L,',(2a)+L;_i(2a) , L^(2a)+L;_, (2a) ^ 



(w±iT)2 + 104o+24-f +f +(2n+l)|eB| (w±ir)2 + 104o+24+f-f+(2n+l)|cB| ^ ' 

x(«4)c 1 _ 1 -ff y/ iy, 4L',_,(2a)+L„(2ff)-L„-i(2a) 4L',_, (2a)+L„(2a)-L„-i (2a) 

•'4 \^) - 2^ '-(„±,r)2 + 104(j+24-f+2nleBl (w±ir)2 + 104o+24+f+2n|cB| J' 

^(^-^V^)-) = y r \\n\ Ll(2a)+Ll^(2a) L;(2a)+£;,_,(2a) 

Jj (CO) 2^ > L(u±!T)2 + 104o+24-f-f+(2«+l)|eB| (oj±iTY + V)Alg+2A\+(+^+(2n+l)\eB\ 

Li(2a)+L;_,(2ff) L,',(2a)+L;_,(2a) , 

(«j±iT)2 + 104^+24 -f+f+(2K+ l)|eB| ~ (w±ir)2 + 104o+24+f-f+(2n+l)|cB| ^ ' 

7(«4).,A ^ 1^-a y / Ll(2a)+Ll^(2a) L;(2a)+Z.;,_,(2a) 

Jfi I"'; 2*^ ^^o*- '-(o)±iT)^ + 10Ajg+2Al-(-^+(2n+l)\eB\ (oj±iT)^+WA]g+2Al+(+^+(2n+l)\eB\ 

Li(2a)+L'__,(2ff) L,',(2a)+L;_,(2a) , 

(bj±ir)- + WAi;g+2Ajf-(+^+(2n+l)\eB\ " (oj±irf + U)Aj,-,+2Ajf+i:-(+(2n+l)\eB\^' 

with being the generalized Laguerre polynomials, a - -k^/\eB\, andx = ^!\eB\. 
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APPENDIX B 

Substituting Eq. (fTTT i into the trace Eq. (14), then it is evaluated after a somewhat tedious calculation 

triys k)yjS k)) = + /r)(w + /r)[Af '(w'Mf V) - Af (w')A^(w) + Af V'Mf V) 

- Af + ieijSgn(eB)(to' ± iT)(co ± ir)[Af'^\co')Af'^'>(co) 

- Af^\oj')Af^\oj) + Af'^\aj')Af'^\co) - Af^\co')Af'^\co)] 
+ {Ikikj - 5ije)[{Bf'^\oj') - Af^\u'))(Bf\oj) - Af^\oj)) 

+ (Bf'^\co') - Af^\aj'))(Bf'^\co) - Af^\aj)) + (fif ^ V) - Af^\c,')) 
X V) - ^ V)) + (Bf V) - Af'\c.'))(Bf'\c.) - Af^\c.))] 

- isgnieB)[6ij(-ir'2kik2 + eij{k] - k]mBf'^\oj') - Af^Hco'))(Bf^\u) - Af ^V)) 
+ <^ V) - Af V) - ^V)) + V) - Af V')) 

X V) - Af ^ V)) + V) - Af V) - ^r'M)] 



/g \^ )h - i£ijSgn{eB)2{ ) [I\ '(w )J\ '( 



where e,j is antisymmetric tensor (en - 1). Integrating over momenta in Eq. (14), we obtain the longitudinal conductivity 

= ^,,(Q 0) = ^^Re da^ ^^^J,^ {X, + - (^)2Z,), (B2) 



where 



(u+ir)2-(10iJo+24)-f-f+|cB| '■i4rw-2f-f+|eB|-2|fB| -2f-f +|eB|-2|fB| J (w+ir)2-(10^|o+24)-^+f+l<'5| 



y r (t^^+r^) (cj+iTf 1 ^ (l-p)(a)-+r ) 

^ L/4rw-2f+f+kB|-2|eB| -2i+^+\eB\-2\eB\i (M+iT)2-(10iJo+24)+f ^ -,4rw-2f-f-|eB|+2|eBI 



-2(-£-\eB\+2\eB\ -iATa>-2(+^-\eB\+2\eB\ -2f +f-|eB|+2|cB| J ' 

_ , . (^+-T)^-(io4o+24)+^ (^+;r)^-(io^lo+24)-^-f+kii| , . . (^+/r)^-(i0i|o+24)-^+f+kfl| 

+ M'4t/'( 2kSj ) + W5l/'( 2kB^ ), 



(B4) 



_ (ai+/r)2-(104o+24)+f+f-kB| (to+/r)2-(104o+24)-f-^+|eij| (o;+/r)2-(10^2o+24)-f+f+kii| 
Zl - ZllA( 2kBi Z2l/'( 2pi ) + 23l/'( ^pj ) 

,, {a,+iTf-(\0Al^+2Al)H-^-\eB\ , {<^+iTf-(\0A\g+2Al)-(-^-\eB\ {a,+iTf-(\0A\g+2Al)+M+\eB\ 
- 2kBi )+Z5l/'( ^j^Bj )-Z6l/'( Jk^j ^ "^^^^ 

,, (^+iTf-(WA\o+2Al)H-^^eB\ (<^+i[f -(mAlo+2A'^)-(^-\eB\ 
+ Z7if/( 2kSj )-Zm 2kfij )' 

with iff{z) being the digamma function and 



2(a>'-+r-) 2(a>+ir)^ (l-l)(^^+r^) (l-l)(a;+/D^ 

-!4rw-2f-f+|cB|-2|cB| -2f-f+|eB|-2|eB| -!4riD-2f+f+|cB|-2|cB| -2f +f +|eB|-2|eB| 

('-pX"-+r-) (i-p)("+'r)- 2(^^+r^) 2(a;+ir)^ 

-!4ra)-2f-f-|cB|+2|cB| -2f-f-|cB|+2|cB| -!4ra)-2f+f-|cB|+2|cB| ~2f+f-|cB|+2|cB| ' 



-!4rw+2f+f-|cB|+2|cB| 2f +^-|cBj+2|cB| ' -i4ra)+2^-^-\eB\+2\eB\ 2f-f-|eB|+21eB| ' 



2(^^+r^) 2(a;+;r)^ _ (i-i)(^^+r^) (i-1)(<^+,t)^ 



^^4 _;4r/,,^9;-_f^l»RI_'>l«RI -);-_^j-UHU9I^HI ' '^S 



/■4ra>+2f-f+|eB|-2|eBl 2f-f +|eB|-2|fB| ' -i4roj+2(+£+\eB\-2\eB\ 2f+f +|eB|-2|eB| ' 
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(w+/r)^-(10^2^+24)+f+f-|cB| («+ir)2-(104o+2-l|)+f+f-l«^l (w+iT)^-(104o+24)+f+f-|<;B| 



-(4r6;-2f-2f+2|cB| 
(w+ir)^-(10i^o+24)+f+f-|cB| 
-2f+2|eB| ' 



-iATb}-2C+2\eB\ 



-2f-2f+2|eB| 



22 



(w+ir)^-(10i2^+2i|)-f-f+|eB| («j+ir)2-(104o+2i|)-f-f+|e£| 



i4rw-2f-2f+2|eB| 
(a)+ir)^-(10i^Q+24)-f-g+|efi| 
-2f+2|eB| ' 



(4rw-2f+2|eB| 



(a;+ir)^-(104o+24)-^-^+kii| 
-2f-2f+2|eB| 



Z3 = 
+ 



_ (bj+iiy -( 1 0/l-, o +2A'-g)-(^+\eB\ 



iAT(ij-2i;+2\eB\ 
(w+ir)2-(10i|o+2i|)-f+f+|eB| 
-2f+2f+2|eB| ' 



(w+iT)2-(104o+24)-f+f+l«^l 
iATo>-l(+24+2\eB\ 



(w+iT)^-(104o+24)-f+f+l«*l 
-2f+2|eB| 



Z4 = 
+ 



(a;+iT)^-(10i^Q+24)+f-g-|efl| 

-(■4rw-2f+2|eB| 
(faj+;T)^-(104o+24)+f-f-|gfl| 
-2f+2f+2|eB| ' 



(w+jT)^ -( 1 04 +24)+f 
-iAToj-l(+24+2\eB\ 



(w+!T)2-(104^+2/l|)+f-f-|eB| 
-2f+2|eB| 



25 



Z6 



(w+/r)--(io,i7,, 



-(4rw+2f+2f+2|cB| 
(&j+tT)^-(104Q+24)-f-f-|eg| 
2f+2|eS| ' 

(0J+iTf-(WAlo+^Al)+CH+\eB\ 

(4rw+2f+2f+2|eB| 
(w+!T)2-(10^2^+2i|)+f+f+|cB| 
2f+2|eB| ' 



(a;+;r)--(iai^^.„+2.ij;)-<.'-f-k-z.'| 

-/4rw+2f+2|eB| 



tt^+/r)--0(),l;;o+2,l„)-^-j;-|gB| 
2f+2f+2|eB| 



(w+!r)2 -( lOiJ^ +2i|)+f +f (w+(T)2 -( 1 Oi^ ^ +2i|)+f 



(4rw+2f+2|eB| 



2f+2f+2|eB| 



Z7 



!4r6j+2<r+2|eB| 
(w+!T)2-( 10^5 o +2^s)+f-f +|eB| 
2f-2f+2|c£| ' 



(w+ir)^-(104o+2^2)_|_^_^^|gg| (w+(T)2-(10iJ,,+2i|)+f-^+|eB| 



;4rw+2f-2f+2|eB| 



2f+2|cB| 



Z8 



(w+(T)2-(104o+24)-f+f-|eB| (w+!r)2-(10^^o+2/l7,)-<r+f-|eB| (w+(T)^-(10/!jo+24)-f+f-|eB| 



-i4rw+2f+2|eB| 
(a;+in^-(104o+24)^^+f-|ei;| 

2f-2^+2|eB| 



-i4rw+2f-2f+2|eB| 



2f+2|cfi| 



The Hall conductivity is then given by 
cr;^(0 ^ 0) 



2c2v||cB|ign(cS) 



'Ay 



(B6) 



with 



X^^ = ai/2-a'i, yH = ^'i/2-(fo'i+fo^), Zh = ci/2-c;, 



(B7) 



where 



a\ = 
+ 



2(l-i) 



2(l + i) 



(«+<T)2-(104o+24)-f+f+kB| -HTuj+2C-^-\eB\+2\eB\ (o,+iTf-(WAj.g+2Al)-C-i+\eB\ -iAToj+2(+i-\eB\+2\eB\ 
u^^t^ r 2(1-7) 2(l+i) 



(«+iT)2-(104o+24)+f [ -(4rw-2f-f-kB|+2kB| ' -i4ru-2(+^-\eB\+2\eB\ 



(B8) 



, (aj+iTf-(lOAlg+2Al)+(. (ui+iVf-(H)A\„+2A\)-i^^+\eB\ (u>MVf-(\0X\g+2A\)-l;,-i+\eB\ 

bi = nnif( 5,;^^ )-n2n )-^lm ) 



2\eB\ ' l^^fy 2\eB\ 

(tj+,r)^-(104o+24)-^+g-|eg| ^ ^ (tj+ir)^-(10ijo+24)-^-f-|eg| ^ 



2kii| 



(B9) 



,, (aj+iTf-(\0A]„+2Al)+{+£-\eB\^ .. (w+TfAmA\f,\-2A\)\C-£-\eB\ (ui^Tf-(WA'-^^,^2A\)-C-£-\eB\ 
C\ = Viif/( )-V2iA( rbiT^ ) + V3(/'( ^ttjt-^^ ) 



(M+iT)--0(),l;„-r2,l„)-^+t-k-Zi| 



) - V5tA(- 



2kB| 

T7(A( §5r^ ) + ysH- 



(ai+iF)- - 


2\cB\ 






2|«-S| , 




-(104o+24)-f+f+le«l 



:|rB| 



2|eB| 



(BIO) 
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= ao + ai M(oJ + /n' - (lOAlo + 2^) - ^ - ^ + \eB\] + a2 MioJ + iTf - (lO^o + 24) - f + ^ + \eB\] 

+ aslnlico + irf -ilOAjg + 2Al) + (], ^ ' 

^, ^ i„r(_ (^^-r)^-(iojo^24)., ^ +^,^(_(^l£f_(||l£4H£) +^3^(_l^l£)!_(g.l^), (B12) 

£>2 - P4lnU 2k3j )+P5fK 2k5j 

„ ,, («+^■^^-(104o+2-^|)-^-f+l«•BI^ , o 1 ((o+ir)2-(l(M2^,+24)-f+f+ki;| 

- Pe'/'l. jklj J+/i7inn jjiSi -> 

+ PS"^*. 2Mi )-P9n 2k5j ) 

, («+ir)^-(104o+24)-f-f-ki}| («+!T)2-(104o+24)-?-f-l««lN 
PlO In 1 (, 2kBi )-Pim 2iiBi 



+ Pi2n 2kBj )-pi3ini{ ; 



/^2 



-i4rw+2f-f-|eB|+2|eB| ' -iATo}+2C+^-\eB\+2\eB\ 



(B13) 



^' _ , p. (a;+ir)^-(104o+24)+f+f-|eB| , (a)+iD^-(104o+2i^)+f-j^-kfi| 

(<u+irf-(lQ4„+24)-f-f-iefii , , ((u+iD^-(104o+24)-f +g-kB| , 

+ r7ini(, j + 7ioini(, j 

_L . ^^rt (w+!T)2-(104o+24)+f+f+kBk _^ , 1„T^^ (o;+,T)^-(104o+24)+f-f+kBk 

+ risinlC jRb^ ) + yi6Jiil( 2k5i ) 

+ ri9inr( ) + r22inr( ) 

(a;+;r)^-(104o+24)+^+f-kfi| ^ , ^_ (a)+ir)^-(104o+2i|)+^-^-k.B| ^ 

+ 72^^ jjiBj J + ren spii 

(a;+in^-(104o+24)-f-g-ki;| , , (a;+,D^-(104„+24)-f+|:-kg| , 

+ rgf^i 2pi ) + ri2(/'( (f^^A^ 

+ risiAi 2Rb^ j + nsiAi 2^5] 

+ r2in 2kB^ ) + y2AH'K 2kBj ; 

+ 2Rg^ ; + r5w '2Ub\ 

, ((^+ir)^-(l(u;o+2,l;,)^f-f-k£| (a)+;T)^-(104o+24)-^+f-kii| 
+ 78<?'(. 2|eB| ;+711V'(. 2\eB\ > 

,, («+ir)2-(io4o+2i|)+f+f+ks| («+!r)2-(io.i^,„+24)+f-f+kB| 

+ y\A9\ 2kSj )+y\i<p\ %B\ 

,, (w+;r)^-(104o+24)-f+f+kB| , («+!r)^-(104^+24)-f-f+kB| 
+ Jn^K 2kBi )+y2a<PK 2kBi 

with 0(z) = In Y{x)dx, and 

/^l ~ -i4rw-2f-f+|eB|-2|eB| -i4rw-2f+f-|eB|+2|eB| -(4rw-2f+f+|eB|-2|eB| -i4rw-2f-f-|eB|+2|eB| ' ^ 

2(l-iK«^+r2) 2(l+i)(«^+r2) 



' F3 _,:4n,i+2,^+/:-l<>BI+2l<.BI ' ^ 



2(i.i)(.^.r^) , 2(i-i)(.^.r^) 

-i4rw+2f-f +|eB|-2|eB| ' ^5 -i4rw+2f +f +|eB|-2|eB| ' ^ ' 



_ T (i^iTf-(WXso+2XKYHH-\eB\ q (cu+ir)^-(10iso+2ifi)'+f+f-kfl| rDio-v 

^1 -i4rw-2f-2f+2|<;B| ^ -i4rw-2f+2|eB| ' ^.^^O-' 

_ ^ (a)+lT)^-(10i5o+2-^«)^+f-f-|eB| o (cu+ir)^-(10A;o+2-<«)^+^-g-|eii| (•9.\C>\ 

^2-^ -iATco-2(+2\eB\ ^ -iATo)-2C+24+2\eB\ ' K^'-^> 

_ ^ (a,+irf-(mAso+2ARf--(-^-\eB\ _ n, (iO+lTf-(WXso+2ARf-(-^-\eB\ CRICW 

■^3 ^ _;4r,,.4.9/4.9^4.9i,>Ri _;4n,.4.9/4.9i»fii ' yDz.\j) 



''4 - ^ -i4rw+2f+2|eB| ^ -Hrw+2f-2^+2|eB| ' >.oz.i; 
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^5 - ^ -i4rw-2f-2f-2|eB| ^ -(4rw-2f-2|eB| ' l,oz,z,; 

_ o (<^+ir)^-(10^so+2^j;)^+f-g+|eji| ^ (<^+ir)^-(10^so+2^«)'+f-g+kg| (TK?^\ 
^6 - ^ -iATo>-2(-2\eB\ ^ -iATb>-2(+^-2\eB\ ' \"^-^) 

_ r, (.a>+iTf-(mso+2ARf-(-i+\eB\ r, (u>+irf-imso+2AR)^-(-i+\eB\ (j,^.. 
n — -iATb>+2(+2^-2\eB\ ^ -iATb>+2(-2\eB\ ' 

_ r, (u>-i-irf-(mso+2ARf-(+i+\eB\ T (<^+irf-ilOAso+2AR)^-(+i+\eB\ m'?'^^ 
^8 ~ ^ -(■4rw+2f-2|eB| ^ -i4rw+2f-2f-2|eB| ' 

2(l+l)(6;+.r)^ (104p+2i|)+^+^-kJ)| 2(l-l)(aj+.r)^ 

^0 - (2f+f-|eB|+2|eB|)2 + 2i+^-\eB\+2\eB\ (,(o+ir)^-aOAl^^+2Ah-(-(+\eB\ (2i-^-\eB\+2\eB\Y 

(mA]o+2Al)H-^-\eB\ 2(l+])(c.+r)^ 

2C-i-\eB\+2\eB\ (^u+iTf -{WAlg+2Al)-(H+\eB\ {2(-^+\eB\-2\eB\f 
(104^+24)-^ + 2(l-i)(..+,rf (104„+24W 

2f-f+kB|-2kB|(a,+,n2-(104o+24)+f (2(+^MeB\-2\eB\f 2f+f+kB|-2kB|(a,+in2_(io4^+24)+f ' 

_ 2(l+l)[(10i|o+2i|)+2f+3f^Sk.B|] _ 2(l-l)[(104„+2^J,)+2f-3f^5k£|] 

^1 - (2f+f-|eJ}|+2|eB|)2 ' ^2 - (2f-f-|eB|+2|ci;|)2 

_ 2(l+i)[(104p+24)-2^-f+kB|-2|efl|] 2(l-l)[(ll)4p+24)-2^+f+|eg|-2[efi|] 
^3 - (2f-f+kB|-2kB|)2 + (2f+f+|cB|-2|cB|)2 

2(l+l)k£|[2(a;+;T)^+2f+f-|g.B|+2|e.B|] 2(l-l)ki}|[2(«;+iD^+2f-^-|eii|+2kg|] 

^'l ~ (2(+^-\eB\+2\eB\f- (2f-f-|eii|+2kB|)2 

2(l-l)|ei}|[2(a,+iD"+2^+f+ka|-2|e£|] 2(l+l)ka|[2(a;+iD^+2f-f+|gg|-2|g.B|] 

+ (2f+f+kB|-2kB|)2 + (2^-f+kB|-2|eB|)2 

(l+IXfc'+'r)^ (l-i;)(io+iTf (X-i^Xaj+iTf (l+|)(a)+iD^ 
"2 - 2(+^-\eB\+2\eB\ 2i-^-\eB\+2\eB\ 2^+f+|eB|-2|eB| 2^-f+|eB|-2|eB|' 



___8kmi+|) 8ki}p(l^i) Sle-Bpd--^) 8kfip(l + ^) 

"3 ~ (2^+f-|eB|+2|eB|)2 (2^-f-|eB|+2|eB|)2 (2f+f+|eB|-2|eB|)2 + (2^-f+|eB|-2|eB|)2 ' 



_ 2a+i^)\eB\[2(cj+irf~2(-^+\eB\-2\eB\] 



_ %(l+'-)\eB\^ (l + i)(«+ir)2 

"5 — r2,>-4./:-l<.BI+2l/.BI12' P6 



(2CH-\eS\+2\eB\y f^" 2(+^-\eB\+2W 



_ 2(\-^)\eB\{2(oj+irf-2(H+\eB\-2\eBW 
~ (2(-i-\eB\+2\eB\f 



8(l-i)kiJp „ (l-i)(a)+(r)^ 



^8 r2V-;-I^RU2I^BIl2> ^! 



(2C-i-\eB\+2\eB\f-' 2^-f-|eB|+2kfl| ' 



_ 2(l-l)|eB|[2(a;+.rp-2f-f-|efl|+2|gfl|] 
"10 ~ (2f+f+|eBh2|eB|)2 ' 

8(l-j)kgP (l-l)(a.+<Tp 
~ (2f+f+|eB|-2|<;£|)2 ' /'IZ ~ 2f+f+|cB|-2|e£|' 



_ 2(l+l)|eB|[2(a;+/rp-2^+f-kii|+2kii|] 
^13 - (2i-^+\eB\-2\eB\f 
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P14 - (2(-^+\eB\-2\eB\)^' "l5 - 2f-f +|eB|-2|eB| ' 



4|e£|[(fai+ir)^-(104o+24)+2^+f-2|efl|] _ 4\eB\[(o)+iTf-(W^lo+2Al)+2(+2^-2\eB\] 
(-2(+2\eB\r i-2C-2i+2\eB\r ' 



_ (u+iT)^-(WAlg+2^l)+(+^-\eB\ (oj+iTf-(m^lg+2Al)+(+^-\eB\ 
~ -2f-2f+2|eB| -2C+2\eB\ ' 



^ = i\eB\' Sk^r 

'3 (-2f+2|eB|)2 (-2(-2^+2\eB\f ' 



4|eg|[(faj+;r)^-(104o+24)+2f-f-2|efl|] _ 4|ea|[(tu+ir)^-(104o+24)+2f-2f-2|efl|] 
(-2f+2|eB|)2 (-2f+2f+2|eB|)2 ' 



'5 (-2f+2|eBD2 (-2f +2f +2|eB|)2 ' 



_ (faj+tT)^-(104o+24)+f-f-|eg| _ (fai+ir)^-(104o+2i|)+f-f-|eg| 
~ -2(+2^+2\eB\ -2f+2|cB| ' 



4|eB|[(«+(T)2-(104o+24)-2f-f-2|eB|] _ 4|eB|[(w+!T)2-(10^2,o+2^|)-2f-2f-2|eB|] 
(2^+2|ei;|)2 (2f+2f+2ki;|)2 ' 



_ 8kfip 8|ei?|- 

/ o (2f +2|cS|)2 (2f+2^+2|eBl)2 ' 



_ (a)+r)^ -( 10--iL +2^^,)-f-g-|ei}| _ (cu+iQ' -( 1 0^»Lj +24)-f-j^-kii| 
y^ ~ 2f+2f+2kB| 2^+2kfi| ' 



4kB|[(a)+(r)2-(104o+2i|)-2^+f-2kB|] _ 4kB|[(«>+ir)2-(104o+2i|)-2f+2^-2kB|] 
(2f+2kfi|P (2f-2#+2kB|)2 



'1' (2f+2kB|)^ (2f-2f +2k£|)^ ' 



_ (w+!rf-(10iJo+24)-f+f-kB| ((j+iT)^-(lOAjg+2Al)-(+^-\eB\ 
y^^ ~ 2f-2^+2k£| 2f+2kB| ' 



4\eB\[(oj+iry-(lQAlg+2Al)+2i:+2f+2\eB\] 4\eB\[(io+iTy-iWAlg+2Al)+2(+^+2\eB\] 
(-2f-2^-2k£|)^ (-2f-2kB|)2 



= 8k-Bp 8|e.Bp 

ru (-2(-24-2\eB\f (-2(-2\eB\f' 

(w+(T)2-(10iJo+2i^)+f+f+kB| (w+(T)^-(104o+24)+f+f+kB| 



ri5 



-2f-2kB| -2f-2^-2kB| 



4kB|[(w+ir)^-(10iJo+2^|)+2f-2f+2kS|] 4kB|[(a)+ir)^-(l(Mjo+24)+2f-f+2kB|] 
(-2f+2f-2kii|)2 (-2f-2kB|)2 
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719 



41eBl[(M+/r)^-(10/l|o+24)-2f-2f+2|cB|] 2|fB|[(w+;r)^-(10/l|o+2/l|)-2f-f+2|e5|] 



720 



8|eBp 



(2f +2^-2|eB|)2 (2(-2\eB\)^ ' 



_ (w+iTf -(104Q+2/l|)-f-f+|eB| (w+irf-(10/l|j,+24)+f-f+|cB| 
- 2(-2\eB\ 2(+2^-2\eB\ ' 



_ 4|cB|[(w+ir)^-(10/!5o+24)-2f+2f+2|eB|] 4|cB|[(w+;r)^-(10/l|o+2/!^)-2f+f+2|eB|] 
'''^Z - (2f-2f-2|eB|)2 (2f-2|eB|)2 ' 



723 



(2f-2f-2|eB|)^ (2f-2|e.B|)- 



_ (w+iT)2-(104o+24)-f+f+kB| (ii>+/r)--(10/l|o+24)-f+f+|cB| 
^24 - 2f-2|cB| 2f-2f-2|cB| " 



The Eqs. ( IB2l i and (IB6l l are further rewritten as Eqs. ( fTSl l and ( fTSI l. 
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